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AN INFRARED DIFFUSE CIRCUMSTELLAR BAND? 

THE UNUSUAL 1.5272 MICRON DIB IN THE RED SQUARE NEBULA 

G. ZasowskF'^ S. Drew ChojnowskF, D. G. Whelan"*, 

A. S. Miroshnichenko’. D. A. Garcia Hernandez®’*^, S. R. MajewskF 


ABSTRACT 

The molecular carriers of the ubiquitous absorption features called the diffuse interstellar bands 
(DIBs) have eluded identification for many decades, in part because of the enormous parameter 
space spanned by the candidates and the limited set of empirical constraints afforded by observations 
in the diffuse interstellar medium. Detection of these features in circumstellar regions, ■where the 
environmental properties are more easily measured, is thus a promising approach to understanding 
the chemical nature of the carriers themselves. Here, using high resolution spectra from the APOGEE 
survey, we present an analysis of the unusually asymmetric 1.5272 ^m DIB feature along the sightline 
to the Red Square Nebula and demonstrate the likely circumstellar origin of about half of the DIB 
absorption in this line of sight. This interpretation is supported both by the velocities of the feature 
components and by the ratio of foreground to total reddening along the line of sight. The Red Square 
Nebula sightline offers the unique opportunity to study the behavior of DIB carriers in a constrained 
environment and thus to shed new light on the carriers themselves. 

Subject headings: circumstellar matter — ISM: lines and bands — ISM: individual objects (Red Square 
Nebula) — infrared: ISM 


1. INTRODUCTION 

For nearly a century, absorption features known as the 
diffuse interstellar bands (DIBs) have been observed su¬ 
perimposed on the spectra of stars, galaxies, and quasars 
along numerous lines of sight through the Milky Way and 
other galaxies (e.g., Heger 1922; Cox 2011). The current 
strongest consensus is that most (if not all) of the distinct 
features are produced by carbon-dominated molecules, 
such as polycyclic aromatic hydrocarbons (PAHs) or 
fullerenes (e.g., Herbig 1995; Sarre 2006), but so far, 
definitive association of any specific species with a spe¬ 
cific DIB has proven elusive. 

The vast majority of, perhaps all, DIB detections have 
occurred in the diffuse interstellar medium (ISM), unas¬ 
sociated with stars (Seab 1995). However, circumstel¬ 
lar environments are particularly interesting potential 
hosts of DIB carriers, because often different types of 
information concerning their temperature, density, and 
chemistry are available than for a particular patch of dif¬ 
fuse interstellar cloud at a generally unknown distance. 
Furthermore, as circumstellar environments are typically 
significantly warmer and denser than diffuse ISM clouds, 
they provide additional, complementary regimes in which 
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to characterize DIB behavior and thus their carriers. 
Searches for DIB features in circumstellar environments 
have been either unsuccessful or inconclusive in a va¬ 
riety of sources, from protostars to planetary nebulae 
(PNe; e.g., Le Bertre & Lequeux 1993; Oudmaijer et al. 
1997; Garcia-Hernandez & Di'az-Luis 2013). Slightly 
more promising candidates have been found by Diaz- 
Luis et al. (2015), who noted absorption features towards 
the PN Tc 1 at wavelengths consistent with DIBs blue- 
shifted within the range of the nebular emission velocity; 
Dahlstrom et al. (2013) and Oka et al. (2013), interpreted 
some unusually broad and asymmetric DIBs towards the 
infrared (IR) source Herschel 36 as most likely due to 
excitation of the DIB carrier molecules by the local radi¬ 
ation field, but see also Bernstein et al. (2015). Broad¬ 
ening and asymmetry may also be caused by the super¬ 
position of DIB features from multiple clouds along the 
line of sight, and many DIBs have been determined to 
have intrinsic substructure and asymmetry (e.g., 6203A 
and I3175A; Porceddu et al. 1991; Rawlings et al. 2014). 

Here we present the discovery of a strongly asymmetric 
instance of the (normally symmetric) 1.5272 ^m IR DIB 
(vacuum wavelength; Geballe et al. 2011; Zasowski et al. 
2015) along the line of sight towards the Red Square Neb¬ 
ula (RSN, also MWG 922; Merrill & Burwell 1949). The 
RSN, a luminous IR source variously classified as a proto¬ 
planetary nebula and a B[e] star (Earners et al. 1998; 
Meixner et al. 1999), has strong hydrogen and metallic 
emission lines that provide constraints on the geometri¬ 
cal and kinematical structure of the nebula (Andrillat 
& Swings 1976; Tuthill & Lloyd 2007; Whelan et al. 
2014). The circumstellar environment hosts a substan¬ 
tial amount of dust, atomic, and molecular material (An¬ 
drillat & Swings 1976; van Diedenhoven et al. 2004), and 
we demonstrate here that it is also responsible for ap¬ 
proximately half of the 1.5272 p,m DIB absorption along 
this line of sight, based on evidence from foreground red- 
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dening and the absorption feature shape and velocity. 

In Section 2, we describe the data used and the extrac¬ 
tion of the RSN’s 1.5272 DIB profile. In Sections 3- 
4, we present the properties of this unusual detection, 
compare it to the typical profiles observed in the ISM, 
briefly examine other IR and optical DIBs along this line 
of sight, and explore plausible interpretations. In Sec¬ 
tion 5, we summarize our conclusions. 

2. DATA 

2.1. H-Band Spectra 

The primary data used are high-resolution, high signal- 
to-noise 77-band spectra of the RSN observed by the 
Apache Point Observatory Galactic Evolution Experi¬ 
ment (APOGEE) survey (Majewski 2012, 2015 in prep). 
APOGEE is a high resolution {R ^ 22,500), 77-band 
(A = 1.51 — 1.70 /rm) spectroscopic survey studying the 
chemodynamical evolution of the Milky Way (MW) with 
a large sample of predominantly red giant stars (Zasowski 
et al. 2013). The survey uses a custom-built 300-fiber 
spectrograph (Wilson et al. 2012) on the 2.5-m Sloan tele¬ 
scope (Gunn et al. 2006); public release of the data began 
in 2013 with the Sloan Digital Sky Survey III (SDSS-III; 
Eisenstein et al. 2011) Data Release 10 (DRIO; Ahn et al. 

2014) and culminated in 2015 with DR12 (Alam et al. 

2015) . 

Because of their usefulness in identifying terrestrial 
atmospheric (telluric) absorption features, several thou¬ 
sand early type stars were also observed by the survey. 
The RSN was observed as part of this sample on 4, 23, 
and 25 May 2013. Figure 1 is a Keck NIRC2 77-band im¬ 
age of the RSN, along with the size and position of the 
APOGEE fiber field of view (red) and some slits (blue 
and orange) used for supplementary information on this 
sightline’s optical DIBs. 


u 

0 ) 

(f) 

o 

o 


o 

(U 

Q 

<1 


10 

5 

0 

-5 

-10 



-10 -5 0 5 

ARA [arcsec] 


10 


Fig. 1.— Keck NIRC2 //-band image of the RSN (grayscale), 
with the apertures of the 2 ” APOGEE fiber (red circle), a 1.2” X 8” 
optical slit (outer blue rectangle), and a 1.0” X 8” optical slit (inner 
orange rectangle). 


All spectra taken with APOGEE are reduced using a 
custom reduction pipeline (Nidever et al. 2015), which 
produces an “apVisit” data file for each stellar observa¬ 
tion. This pipeline calculates the radial velocity (RV) 
for each visit, but as it is not optimized for emission 
line spectra, we calculate our own RVs with a more ro¬ 
bust treatment of the spectral features and use those to 
stack the apVisit data into a single source spectrum (with 
S/N > 460). 

2.2. Extraction of 1.5272 pirn DIB Profile 

The 1.5272 fim DIB overlaps the Brackett-19 hydro¬ 
gen emission line (Brl9; Aq = 15 264.725 A). Because 
the DIB appears very close to where the Brl9 line drops 
to the continuum, careful treatment of the intrinsic Brl9 
profile shape is required. Here, we model the Brl9 pro¬ 
file empirically using all of the unblended hydrogen lines 
available in the APOGEE spectrum (i.e., Bril, Brl3, 
Brl6, Brl7, and Br20). In the rest velocity frame, we 
calculate a quadratic relationship between the observed 
flux F\{Av) and the Aq at each Av of the Brackett lines. 
Then, we interpolate this relationship at the Brl9 rest 
wavelength to construct our empirical Brl9 profile. Fig¬ 
ure 2a shows these profiles (colored lines) and the inter¬ 
polated BiT9 model (black dashed line), along with the 
observed BiT9-|-DIB profile (thick blue line). We mea¬ 
sure a heliocentric RV of 37.8 km s“^ from this line, cor¬ 
responding to a velocity relative to the Local Standard 
of Rest (ulsr) of 49.8 km s“^. 

3. ANALYSIS 

Figure 2b shows the strong residual absorption {R\) 
remaining after normalization by this empirical Brl9 pro¬ 
file model (with two fitting scenarios that are discussed 
in Section 3.2). To our knowledge, there are no other 
atomic or molecular transitions at this wavelength that 
could be strong and broad enough in this source to pro¬ 
duce the absorption observed. 

We calculate an integrated equivalent width Wdib = 
557 mA by summing (1 — R\) in the window —150 < 
flsr a 250 km s“^; random uncertainties based on the 
“flat” span of R\ between —300 < ulsr < —150 km s“^ 
are <1 mA. Assuming the mean ratio Wdyb/Av = 
102 mA mag ^ for this DIB (Zasowski et al. 2015), 
and a typical Rv = 3.1 extinction law, this equivalent 
width corresponds to a line-of-sight dust reddening of 
E{B — V) « 1.75. At high reddening, the dispersion in 
the Wdib/Av ratio is <40 mA mag“^. 

Of course, for any circumstellar material, the extinc¬ 
tion law (and Wbib/Av ratio) may be atypical. How¬ 
ever, we can make some independent estimates of the 
interstellar reddening. The optical DIBs also measured 
as part of this analysis suggest a total reddening of 
E{B — V) ^ 1.0 — 1.8 mag, with a dispersion of about 
0.5 mag, depending on the optical spectrum used (see 
Section 3.4). The foreground-only reddening, as mea¬ 
sured from PanSTARRS-|-2MASS photometry by Green 
et al. (2014, 2015) for the T x T region around the RSN’s 
position, is E{B — V) « 0.57llo;o| and « 0.87llo;o5 mag 
at heliocentric distances of ~2 and ~3 kpc, respectively. 
These lower values are consistent with the presence of 
additional dust along the line of sight being traced by 
the IR DIB, and perhaps some of the optical ones. 
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Fig. 2.— Measurement of the overlappping Brl9 and 1.5272 fim DIB profiles, (a): The unblended Brackett profiles (thin colored lines) 
are used to interpolate the intrinsic Brl9 profile (dashed black line). The observed Brl9+DIB profile is shown in blue, (b): The residual 
DIB absorption, versus velocity ulsri with single (top, red) and dual (bottom, green) Gaussian fits. The vertical lines indicate the centers 
of the Gaussian components, and the gray arrow is the i>lsr of tho RSN’s Brl9 line. 


3.1. Asymmetry of Profile 

To quantitatively compare the asymmetry in the RSN’s 
1.5272 fj.m DIB to that of the typical 1.5272 /tm DIB fea¬ 
tures observed in the ISM, we use a simplified version of 
the line bisector estimator (e.g., Toner & Gray 1988), 
which compares the midpoints of a profile near its top 
and bottom (15% and 85% of the amplitude; left panel of 
Figure 3). A difference in these midpoints (here, 6v) in¬ 
dicates asymmetry in the prohle. When compared using 
this metric to a set of several thousand 1.5272 /im DIB 
features that sample the ISM towards cool giant stars® at 
identical spectral resolution, the RSN sightline appears 
more asymmetric than all but ~2% of the larger sam¬ 
ple (right panel of Figure 3). Note that this comparison 
set does contain noisy sources, because restricting it to 
only those sightlines with “good” DIB fits would explic¬ 
itly remove any asymmetric features. Visual inspection 
of those ^2% asymmetric sources confirms that the vast 
majority are contaminated by residual stellar features, 
though we cannot rule out a few potentially intrinsically 
lopsided DIBs that will be good candidates for follow¬ 
up observations. None of the sightlines with potentially 
asymmetric 1.5272 /im DIB features lie closer than 2.4° 
to the RSN, and none of the DIBs in the ~60 sightlines in 
the dataset at closer projected distances show any signs 
of asymmetry. 

3.2. DIB Fit Components 

® This comparison set comprises all sightlines from Zasowski 
et al. (2015) with cr(Rx) < 0.03, A^visits = 1 or Wgcatter < 
0.5 km s“^, a{Rx)/rj{Fx) < 0.55, and A^pix > 1350; see Section 2.2 
of that paper for details. 


We turn now to different options for fitting this fea¬ 
ture. The upper Rx profile in Figure 2b is overplotted 
in red with the Gaussian that best reproduces the blue 
(lower velocity) side of the feature; the extended red ex¬ 
cess is obvious. (Fitting the entire feature with a single 
Gaussian yields significantly noisier residuals.) The bulk 
of the feature is well-represented by this Gaussian with a 
FWHM of 3.5 A. If this DIB feature has a single origin, 
the red excess may be due to changes in the DIB car¬ 
rier’s electronic properties (Section 4). The fact that the 
deepest part of the feature appears at ulsr ~ 12 km s“^, 
nearly 37 km s“^ offset from the nebular Brackett emis¬ 
sion, is a strike against its origin in the nebula. The 
mismatch does not completely rule out a circumstellar 
origin, as a range of velocities (of tens of km s“^) is ob¬ 
served in other absorption and emission lines definitively 
associated with the nebula, arising from the superposi¬ 
tion of multiple nebular components (e.g., Whelan et al. 
2014). However, a simpler solution may be sought. 

The preferred scenario is shown with the lower Rx pro¬ 
file of Figure 2b, overplotted in green with the best simul¬ 
taneous fit dual-Gaussian components, which together 
reproduce the entire feature very well. Each component 
corresponds to an approximately equal Wdib contribu¬ 
tion of 277 mA, and thus an E{B — V) « 0.87 mag. 
This is consistent with the foreground reddening pre¬ 
dicted by the PanSTARRS-|-2MASS map in the distance 
range of the RSN, implying an additional reddening of 
0.87 mag due to circumstellar material. Using instead 
a dense environment extinction curve of Ry = 5 with 
one component, this additional reddening is closer to 
E{B — U) ~ 0.54 mag, which brings the total redden¬ 
ing to 1.4 mag. Both of these are within the (albeit 
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Fig. 3.— Left: The measurement of the asymmetry estimate Sv for the RSN’s DIB. Right: The distribution of ISM Sv values, compared 
to that from the RSN sightline. 


broad) range of reddening spanned by the optical DIBs, 
suggesting that some of the optical DIB absorption may 
arise from the circumstellar medium as well. 

The use of reddening arguments alone as evidence 
for circumstellar DIBs is notoriously tricky, because of 
the large uncertainties and intrinsic dispersions in the 
DIB/£’(i? — V) relations (e.g., Luna et al. 2008). Thus, 
further critical evidence for this scenario lies in the line- 
of-sight velocities of the two components. The narrower 
feature (FWHM = 3.1 A) has ulsr = 10.3 km s“^, 
well within the range spanned by the optical DIBs and 
the interstellar Na DI and K I lines. In contrast, the 
wider feature has rinSR = 49.6 km s“^, nearly identical 
to that measured for the Br emission lines of the RSN 
itself. Under the interpretation of a circumstellar ori¬ 
gin, the larger width (FWHM = 8.7 A) may be seen as 
simply velocity broadening due to turbulence or rotation 
within the nebula. If the foreground component approxi¬ 
mately represents the intrinsic feature width, this implies 
^47 km s“^ of broadening within the nebula, consistent 
with the range of velocities observed in different nebular 
components (e.g., Whelan et al. 2014). Thus, if the two 
postulated components indeed represent different phys¬ 
ical origins, then they serve as strong evidence for DIB 
carriers in the circumstellar environment of the RSN. 

3.3. Other H-Band DIBs 

Ten known DIBs lie within APOGEE’S IR wavelength 
coverage (Geballe et al. 2011; Gox et al. 2014), includ¬ 
ing the 1.5272 /im feature. Based on the small num¬ 
ber of sightlines in which the others have been observed, 
their strengths are expected to be at most ~30% of the 

1.5272 /im DIB, if detected at all. Despite the strong 

1.5272 /im absorption in the RSN, however, the diffi¬ 
culty of measuring these weaker bands is greatly exac¬ 
erbated (indeed, dominated) by the complex spectrum 
of the RSN — the numerous strong, blended, and some¬ 
times unidentified emission lines make determination of 
the underlying continuum extraordinarily tricky. (The 

1.5272 /im DIB lies in an exceptionally “smooth” region.) 
Figure 4 demonstrates this with six of the DIBs. These 
features are indeed quite weak, but we do not see strong 
evidence for any feature asymmetries. 

3.4. Optical DIBs 


We also studied two high-resolution optical spectra, 
to assess the properties of optical DIBs and interstellar 
lines in the context of our proposed scenario. One spec¬ 
trum (hereafter, “MDO”) was observed with the Har¬ 
lan J. Smith 2.7 m telescope at McDonald Observatory 
on 6/7 September 2009, using the 1.2" slit shown in blue 
in Figure 1. The R = 60,000 spectrum, which spans 
4787-8861 A and has been smoothed with a moving box¬ 
car filter (with a kernel of 11 points), has a S/N of ^20 
near 4860 A and ~90 near 6563 A. The other (hereafter, 
“UVES”) was observed with the ESO VLT-I-UVES on 
27 August 2009, using the 1" slit shown in orange in 
Figure 1. The UVES spectrum spans 5655-9464A with 
R = 42, 000 and a median S/N of 47.2. Of the DIBs that 
are in regions relatively free of contaminating lines, we 
could reliably fit seven features using a simple Gaussian, 
shown in Figure 5 (MDO on the left, UVES on the right). 

In comparing the two spectra, the largest differences 
are in the NaD doublet, which is saturated in the UVES 
spectrum but not in the MDO one, and in the strength 
of the 5780A and 5797A DIB features (bottom curves 
in Figure 5). The UVES 5780A and 5797A DIB equiv¬ 
alent widths are approximately twice as high as their 
counterparts in the MDO spectrum, and the other DIB 
equivalent widths are between 1.3-1.9x as high. The 
slit widths used differ only slightly (1.2" for MDO and 
1.0" for UVES), and the slits were centered on the same 
position to within the uncertainty of the telescope point¬ 
ing reported in the file headers. To assess the possi¬ 
ble impact of small angular offsets sampling small-scale 
variations in the ISM across the nebula, we also visu¬ 
ally examined the NaD and K I lines from three optical 
VLT-I-X-Shooter spectra^° separated from each other by 
1 — 2" on the sky. The spectral resolution is lower, but 
we see no difference in the shape or strength of these in¬ 
terstellar features. More detailed mapping of any ISM 
structure or DIB variability across the face of the nebula 
is reserved for future work; here, we take the consistency 
of the X-Shooter spectra as evidence that the discrepan¬ 
cies between the MDO and UVES spectra are not due 
to strong variations in the ISM density on the scale of 

These data, spanning 534-1020A, were obtained from the ESO 
archive. One spectrum were observed on 02 July 2010 with median 
S/N = 123, and the other two on 07 July 2012, with median S/N = 
145 and 144. 
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Fig. 4.— Examples of other DIBs in APOGEE’S wavelength range. Top: Observed flux (black) and polynomial continuum (green). 
Bottom: Continuum normalized flux, with DIB central wavelengths and FWHMs (black brackets; average FWHMs from Geballe et al. 
2011; Cox et al. 2014). Pixels flagged as BAD or with high errors in any of the apVisit flux arrays have been masked. 


arcseconds. For the remainder of this analysis, we will 
quote values derived from the two spectra separately. 

Qualitatively, the spectra are very similar. The criti¬ 
cal point is that none of the DIB features, in either opti¬ 
cal spectrum, display the strong asymmetry seen in the 
1.5272 /im DIB (except possibly 6203A, which has been 
previously noted as intrinsically asymmetric; Porceddu 
et al. 1991). We attempted to fit two simultaneous Gaus- 
sians to these features, emulating the dual components 
seen in the 1.5272 p,m DIB, but the resulting fits were 
quantitatively no better than those with a single Gaus¬ 
sian (including for 6203A). We note that at these wave¬ 
lengths, the shift expected for a 37.8 km s“^ offset is 
^0.7A, which is discernible even in these broad (and 
potentially substructured) features; however, there are 
degeneracies among the fit parameters (e.g., amplitude 
ratio vs. velocity broadening), and thus we cannot con¬ 
clusively rule out the presence of red excess absorptions 
with strength ratios or velocity properties slightly differ¬ 
ent from those of the IR DIB components. 

The total equivalent widths of the optical features 
(Wopt), combined with the Wopt/Wi.5272pm relations cal¬ 
culated by Cox et al. (2014), predict an absorption of the 
1.5272 /im DIB that is ~30 — 50 ± 10% what is observed 
(from the MDO and UVES data, respectively); we note 
that these relations are based on a small number of sight¬ 
lines, which probe interstellar material with potentially 
different chemistry or radiation exposure than that to¬ 
wards or in the RSN. 

Using the mean Wopt/E{B — V) relations from Jen- 
niskens & Desert (1994), we calculate a reddening of 
E{B - U) ~ 1.0 ± 0.2 mag (MDO) and E{B - U) ~ 
1.8 ± 0.5 mag (UVES). These are larger than the fore¬ 
ground reddening from the PanSTARRS-|-2MASS map 
at 2-3 kpc, approaching that of the total estimated from 
the 1.5272 /im DIB. This is consistent with a circum- 


stellar contribution to the DIB absorption, on the or¬ 
der of 15-50% of the total. This result on its own is 
merely suggestive, and we emphasize the large dispersion 
in the W^pt/ E{B — V) relations (both intrinsic scatter 
and that due to measurement uncertainties), as noted in 
Section 3.2. However, when combined with the velocity 
comparisons below, it is worth noting the consistency of 
the reddening with excess DIB absorption as measured 
in several DIBs. 

Adopting the central wavelength Ac values of Hobbs 
et al. (2008), we find DIB ulsr values spanning 2 — 
20 km s“^, centered on the value measured for the “fore¬ 
ground” 1.5272 /im DIB (wlsr = 10.3 km s“^), the 
interstellar 5895 A Na D1 line in the MDO spectrum 
(flsr = 8.5 km s“^), and the interstellar 7701 A K I line 
in the UVES spectrum (<ulsr>= 15.2 km s“^). The 
optical DIB feature widths are also consistent with the 
widths from the interstellar cloud studies of Jenniskens 
& Desert (1994) and Hobbs et al. (2008), with some scat¬ 
ter of <8%; as noted above, though a dual-component fit 
does not produce a quantitatively better match to these 
features, we cannot rule out the possibility that a small 
fraction of circumstellar material may be contributing to 
the observed widths. 

4. OTHER INTERPRETATIONS 

There exist other potential explanations for the “addi¬ 
tional” long-wavelength 1.5272 /rm DIB absorption, but 
all face serious problems. One interpretation is that the 
ISM’s velocity distribution towards the RSN has a pecu¬ 
liar structure, such that the “high” ulsr feature is due to 
an unassociated ISM cloud. In a two-cloud scenario, this 
would require the second one to have nearly precisely the 
same radial velocity as the RSN and an internal veloc¬ 
ity dispersion of >47 km s“^, neither of which is readily 
apparent in the observed profiles of the optical DIBs. 
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Fig. 5.— Single Gaussian fits (blue lines) to seven optical DIBs. The normalized spectra are offset by increments of 0.5 for visual clarity. 
Left: Smoothed R = 60, 000 spectrum from McDonald Observatory. Right: Original (gray) and smoothed (black) R = 42, 000 spectrum 
from VLT+UVES. 


To assess this in greater detail, we examined the ve¬ 
locity structure of the 7701 A K I line from the UVES 
spectrum (Section 3.4). This is shown in the top panel 
of Figure 6, along with the best three-Gaussian fit rep¬ 
resenting the three resolved subpeaks. In the bottom 
panel, we “model” the 1.5272 /rm DIB with these three 
components, fixing the central velocities, velocity disper¬ 
sions, and relative amplitudes, but scaling the shortest- 
wavelength amplitude peak to match the deepest DIB 
absorption. These components (and their summed to¬ 
tal) are shown in red in the bottom panel of Figure 6, 
compared with the dual-component best fit presented in 
Figure 2b and Section 3.2. It is clear that while the 
K I profile does reflect ISM components at higher line- 
of-sight velocities, these are not sufficient to explain the 
long-wavelength 1.5272 p,m DIB absorption, under the 
assumption that the DIB carrier and K I are mixed in 
the ISM (like carriers of other DIBs are; e.g, Galazutdi- 
nov 2005). 

Another possible explanation is the presence of an 
additional, hitherto undetected DIB blended with the 
1.5272 fim feature (e.g., as postulated by Bernstein et al. 
2015, for the Her 36 asymmetric DIBs), or intrinsic asym¬ 
metry in this feature. However, this is not supported by 
any other observation of this feature to date, including in 
heavily reddened, early-type stars with strong 1.5272 fim 
features that are symmetric at a high level of confidence. 


Some examples are shown in the right-hand panels of 
Figure 7 (7b and 7d). These spectra are also APOGEE 
data and have a resolution identical to that shown in 
Figure 2. The red lines are single Gaussian fits. If the 
long wavelength extension of the RSN’s 1.5272 /im DIB 
were a result of an intrinsic asymmetry or a distinct DIB 
feature, it should be apparent in these spectra as well, 
but all sightlines show a clear symmetry. 

We confirmed that the asymmetry is not an arti¬ 
fact intrinsic to the APOGEE data by comparing the 
Brl9-|-DIB profile to that in an X-Shooter spectrum^ 
we find that they share the same shape. Figure 7 in full 
presents two additional examples of extracting 1.5272 fim 
DIB features from the APOGEE spectra of stars with 
Brl9 emission, using the methodology described in Sec¬ 
tion 2.2. In the same style as Figure 2a, Figure 7a shows 
the weak, double-peaked Br emission profiles of ABE- 
022 (Chojnowski et al. 2015), with the DIB clearly visi¬ 
ble. Figure 7c shows the same for ABE-137, a star with 
single-peaked emission profiles more like the RSN. The 
right-hand panel of each row contains the residual ab¬ 
sorption after removal of the interpolated BiT9 profile as 
black lines, along with a single Gaussian fit (red lines). 
For comparison, we also reproduce the RSN’s 1.5272 fim 

Obtained from the ESO archive, spanning 994-2479A, ob¬ 
served on 02 July 2010 with a median S/N = 538. 
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Fig. 6. — Extrapolating a multi-component 1.5272 pm DIB from 
the K I ISM line. Top: The UVES spectrum towards the RSN, 
showing the clearly structured K I line at 770lA. The blue lines rep¬ 
resent the best-fitting Gaussian components (dotted) and summed 
total (solid). Bottom: The observed 1.5272 pm DIB (black), with 
the best-fit components presented in Figure 2b and Section 3.2 
(green). The red dashed lines are simulated DIB components with 
the same central velocities, velocity dispersions, and relative am¬ 
plitudes as the K I components above. The solid red line is their 
sum. 

DIB in Figure 7b in blue, and in Figure 7d, we over- 
plot in green the normalized spectrum of ABE-Q09 (also 
Chojnowski et al. 2015), a very early type star with 
no emission lines to complicate the DIB measurement. 
The isolated DIB absorptions are highly symmetric, with 
shapes that are nearly perfectly Gaussian and that are 
extremely similar to the DIB profile in the emission-free 
spectrum. This demonstrates clearly that the asymme¬ 
try in the RSN’s 1.5272 /xm DIB is not an artifact of our 
deblending process. 

We cannot fully rule out the possibility that the RSN’s 
1.5272 fim DIB feature arises entirely from circumstellar 
carriers strongly impacted by their local environment. 
Dahlstrom et al. (2013) reported the discovery of optical 
DIBs with pronounced red wings along sightlines towards 
Herschel 36, a system of O stars illuminating the MS HII 
region and in close proximity to the Her 36 SE strong 
IR source. Qualitatively, these wings are very similar to 
the one observed here in the RSN (see, e.g., their Fig¬ 
ure 6). Chemical models calculated by Oka et al. (2013) 
demonstrate that these wings can be produced when a 
strong radiation held drives some electrons in the DIB 
carrier molecules into excited states. However, the need 
to invoke an outhow of nearly exactly the same veloc- 
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Fig. 7.— Extraction of the 1.5272 fim DIB feature from two 
additional emission-line stars (ABE-022 and ABE-137; Chojnowski 
et al. 2015). Panels (a) and (c) are identical in form to Figure 2a. 
Panels (6) and (d), like Figure 2b, show the residual absorption 
in these spectra after the Brl9 feature has been removed (black 
lines), along with a single Gaussian fit (red lines). In panel (6), 
we also show the RSN’s DIB profile (blue line) for comparison. In 
panel (d), we scale and overplot the spectrum of the early type star 
ABE-Q09 (green line Chojnowski et al. 2015). See text for details. 

ity as the RSN’s motion relative to the Sun, such that 
the main feature is centered at the same velocity as the 
foreground ISM, renders this a less likely interpretation. 

5. SUMMARY AND CONCLUSIONS 

We have presented the detection of a peculiar, highly 
asymmetric instance of the 1.5272 /xm DIB along the 
line of sight towards the Red Square Nebula, using high 
resolution 7J-band spectra from the APOGEE survey. 
Strongly asymmetric profiles have never yet been ob¬ 
served in this particular DIB feature, or any at longer 
wavelengths. 

The likeliest explanation is that the feature towards 
the RSN is a blend of absorption arising from the fore¬ 
ground ISM and from the nebular material itself. This 
interpretation is supported by the radial velocities of the 
two blended components, which align nearly perfectly 
with the foreground ISM (e.g., as traced by interstellar 
Na or K) and with the RSN’s hydrogen emission lines, 
respectively. In addition, the amount of interstellar red¬ 
dening “predicted” by the interstellar DIB component 
nearly exactly matches the foreground reddening traced 
by high-quality three-dimensional dust maps. The total 
reddening predicted by the whole IR DIB feature (in¬ 
terstellar plus circumstellar) is roughly comparable than 
the reddening predicted by the optical DIBs, suggest¬ 
ing that at least some of the optical DIB carriers re¬ 
side in the nebula itself or that the relationship between 
the 1.5272 /xm DIB and the optical DIBs, currently only 
coarsely constrained, is substantially different in the cir¬ 
cumstellar material. 

The asymmetric 1.5272 /xm absorption feature ob¬ 
served towards the Red Square Nebula is likely the most 
conhdent detection to date of a diffuse circumstellar band 
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(DCB; see also Diaz-Luis et al. 2015), and if genuine, is 
also the first detection of an IR DCB, cleanly separated 
from the ISM along the same line of sight in both ve¬ 
locity and dust reddening. Because the circumstellar en¬ 
vironment (its radiative properties, temperature profile, 
geometry, etc.) can be understood and disentangled, this 
source will become a powerful laboratory for constraining 
the response of the DIB carriers to these environmental 
factors, and thus constraining the nature of the carriers 
themselves. 
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